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Abstract- This paper is concerned with change point detection 
using recursive estimation in non-linear time series. Specially, we 
consider the Black and Scholes (1973) option pricing model 
under stochastic volatility assumption. The test statistic is given 
and we describe how its null distribution may be simulated. 
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I. 



INTRODUCTION. 



A key assumption for making inference using a statistical 
model is the model stability over time. This condition, 
however, may be violated in practice and all statistical outputs 
become wrong. For example, in financial time series setting, 
the parameters of a GARCH series may be shifted, because of 
changes in economic environment and forecast results fail. In 
this case, a structural break happened and a common 
technique to assess the constancy of a model is change point 
analysis. Change point detection in financial time series is 
very important, in practice. For example, the portfolios 
volatility may increase as risk premium rises and a change 
point analysis is needed. Kim et al. (2000) considered the 
problem of multiple change points in GARCH models. 
Hillebrand and Schnabl (2003) studied change point detection 
in volatility of Japanese foreign exchange rates under 
GARCH modeling. Lee et al. (2006) detected change points 
in a diffusion process using one step estimators. Iacus (2008) 
considered the change point problem in volatility of a 
diffusion process. An excellent reference in change point 
problems is Csorgo and Horvath (1997). 

Both recursive estimation and recursive errors topics play 
important roles in statistical inference, especially in change 
point analysis. Plackett (1990) and Bartlett (1991) studied the 
recursive least square estimations in a linear regression model. 
Hannan (1980) considered recursive estimators based on 
ARMA models, for more details see Hannan (1980) and 
references therein. The logic behind the recursive estimation 
is to compute parameter estimates with a fraction of sample. 
If the actual parameters are fixed over the time, then these 
estimates are very close together. If a change (or some) has 
happened, then the recursive estimates will show this shift. 
Brown et al. (1975) used an alternative method. They derived 
estimated recursive residuals in a regression model and 
applied the cusum method. Their approach is used by Kramer 
et al. (1988) in dynamic models and it is continued up to this 

time. Throughout the current paper, W t and B t are two 
standard Brownian motions defined on te[0,T]such that 
cor(dw t ,dB t ) = p 

Hannan (1980) presented a regression model 
representation of ARMA process and used Plackett algorithm 
to derive the recursive estimates of regression coefficient. As 
follows, we review this method with two new models. 



Example 1. Let X t be a ARCH(q) model, i.e., X t = \\S t 
where S t ' S are independent and identically distributed 
variables with Efa) = and E(s 1 ) = 1 and 

Let rj t = S t — 1. Then following Chatterjee and Bose 



(2005), we have y t = Z ( p + h^7] t 



where 



/3 = {/3 Q ,/3 1 ,...,P q y and z;={l,Xl v ...,Xl q ) and 

y t = X t . Therefore, following Brown et al. (1975), the 
weighted recursive estimates of jB is 

br=br-i + A"WCy,-2,Vi)A = o, 

And A,, is the matrix whose i-th row is z. satisfying in 
the following relation 

A-l _ A-l i K A-l Z r Z r A--1 

i + h; 2 z r A;\z r 

Also, the cusum test procedure of recursive estimated 
residuals e r may be applied to detect a change point, where 

Recursive estimation also appears in maximum likelihood 
estimation. However, in the next example, we study it for 
diffusion processes. 

Example 2. Let X t be a continuous time AR(1) model 
given by 

dX t =aX t dt + adX t . 
Using the Euler method for a small bandwidth h , we 
have Y. = aX. + a*Z. , where 

h H '/ 

And Z t are independent and identically distributed come 

form N(0, 1) distribution. Again, the Brown's method may 

be applied here to detect change points. In the next section, 
we see how the recursive estimation may be applied for 
Black-Scholes formulae. 



IJED Vol. 2 No. 1, 2012 PP. 13-15 © 2011-2012 World Academic Publishing 



13 



International Journal of E-Business Development (IJED) 



II. BLACK-SCHOLES MODEL 

Recursive estimation also appears in continuous time 
series. For example, Lee et al. (2006) proved that the limiting 
distribution of the one-step ML estimator converges to 
Brownian motion as sample size goes to infinity. Kutoyants 
(2004) derived an recursive relation for ML estimator in 
diffusion processes. This estimation procedure plays a central 
role in many fields of applications such as economic, 
engineering and finance (see Ljung et al. (1988) and Young 
(2011)). Test statistics in change point field (say cusum 
statistic) are constructed based on the RML (see Lee et al. 
(2006)). One of useful financial models is Black and 
Scholes's (1973) option pricing formula. This model has had 
a huge influence on growth of financial mathematics, say, for 
making decision in hedging or valuing European call and put 
options. This formula implies that the percentage change in 

price of a stock S t during life time dt (too small) follows a 

normal distribution. More precisely, it states that 

ds t = <ps t dt + <j t s t dw r 

Here, we assume that ^C is an unknown parameter which 
should be estimated. There are many estimation methods for 
cp such as the recursive maximum likelihood (RML) method. 

In fact, it is the maximum likelihood estimation (f> t based on 

observations x,se[0,t]. The standard Black-Scholes 

model assumes that the underlying volatility is constant over 
the life of stock. However, this assumption isn't true, in 
practice. The stochastic volatility models are solutions to this 
difficulty. These models consider the underlying stock's 
volatility as a random process, i.e., 

da\ = juafdt + £crfdB t . 
In this case, it is easy to see that 

y t - ( p = C- 1 ^cj- s l dw s , 

f T -2 

Where C T = CT du . Therefore, the conditional 

1 Jo " 

distribution of (f> t —(p given C T is normal distribution 
N(0,C?). 

In this paper, we test the existence of a change point in 
parameter (p based on RML estimator (f> t . This means 
parameter (p is changed at some unknown time points 
t from <p x (for t<t Q ) to <p 2 (for t>t Q ). To this end, 

following Kutoyants (2004), the maximum of likelihood 
function is the maximum of following process with respect to 
argument t 



U = min U T (t) 



o<t<r 



r 

Jo 



(ds u -p lt s u duf rT(ds u -p 2t s u du) 



- + 
cr'.s; •" a's 



rV 

u u 



J 



V 

u u 



Here (p lt and (p 2t based on observations {x s } SE r ti anc ' 
{x s } s r tr , , respectively. The null hypothesis of no change 

point is rejected when the above test statistic is large. 
Following Bai (1994), we can show that the test statistic is 
reduced to 



= min C T a c (1 - a t ) |(^ lt - p 2t f , 

-1 f l -2 

Where CL=C T \(J du . It can be shown that, under 

i J Jo " 

the null hypothesis, the test statistic is minimum of absolute 
value of the following process 

* [X t -a t X T l 

y]C T a t Q.-a t ) 

With X, = { cr _1 dW . The X-d f X T is the diffusion 

[ Jo 

type of X t . Therefore, we can simulate the null distribution 
of test statistic by simulation of Bessel-type tied-down 

process of X. = \ (T~ dW„ . 

Jo 

Remark 1. Under the alternative hypothesis, there is a 
change point and (p x ^ (p 2 . It is easy to see that 



\(Vu 'Pit) - £ (^ir -&2t)\ ec l uals t0 



This shows that (cp u — (p 2t ) is very close to its 
expectation. Therefore, it behaves similar to E{ff> xt —<p 2t ) ■ 
Let the magnitude of change be 8 = <p x — <p 2 . Without loss 
of generality, suppose that 8 is positive. One can check that 
the conditional expectation of E((p u —(p 2[ ) given 



cT [0ir] =(cy t ,0<t<T) 



l-a t 
1-a, 



for t < t„ and 



for t>t n . Therefore, 



-E(U T (0) 



Ei^d-ajl^-) 



Ei^aJ^-) 



t<t n 



t<t . 



The function at the right side of above equation is 
increasing before t and is decreasing after t This shows 

that the minimizing of U T (t) is a consistent estimator for 

change point t Q . 

III. CONCLUSION 

In this paper, we propose the use of recursive estimators 
for change point detection. Our approach is constructed by 
representing the non-linear models as regression type model. 
We consider this estimator for GARCH and diffusion series. 
As an important special case, the Black-Scholes model is 
considered and the null distribution of test statistic is 
simulated. This paper can be considered as generalization of 
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recursive analysis in discrete time-series to continuous time 
stochastic processes. 
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